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Abstract. The photoinduced 3D orientational structures in films of a liquid crystalline polyester, containing 
azobenzene side groups, are studied both experimentally and theoretically. By using the null ellipsometry 
and the UV/Vis absorption spectroscopy, the preferential in- plane alignment of the azobenzene fragments 
and in-plane reorientation under irradiation with linearly polarized UV light are established. The uniaxial 
and biaxial orientational order of the azobenzene chromophores are detected. The biaxiality is observed in 
the intermediate stages of irradiation, whereas the uniaxial structure is maintained in the photosaturated 
state of the photoorientation process. The components of the order parameter tensor of the azobenzene 
fragments are estimated for the initial state and after different doses of irradiation. The proposed theory 
takes into account biaxiality of the induced structures. Numerical analysis of the master equations for the 
order parameter tensor is found to yield the results that are in good agreement with the experimental 
dependencies of the order parameter components on the illumination time. 
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1 Introduction 

Effect of photoinduced anisotropy (POA) implies that the 
optical anisotropy revealed itself as dichroism of absorp- 
tion or birefringence is brought about in medium under 
the action of light. The capability of having the light- 
controlled anisotropy makes the materials that exhibit 
POA very promising and highly perspective for use in 
many photonic applications such as optical data storage 
and processing, telecommunication and reversible holog- 
raphy [i). In addition, it was found that substances with 
POA effect serve as excellent aligning substrates for liquid 
crystals |, I . 

Polymers containing covalently linked photochromic 
moieties such as azobenzene derivatives are known as azo- 
polymers. These materials exhibit POA of extremely high 
efficiency: the value of photoinduced birefringence in azo- 
polymers can be as high as 0.3 and the dichroic ratio of the 
absorption is over 10. It makes azopolymers particularly 
suitable for the investigation of light induced ordering pro- 
cesses. This is why in the last decade these polymers have 



been the subject of intense experimental and theoretical 

studies §, § 0, r, @, a ecl nil m. 
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The accepted mechanism of POA induced by the lin- 
early polarized UV light involves induced trans-cis-photo- 
isomcrization and subsequent thermal and/or photochem- 
ical czs-^rans-back-isomerization of the azobenzene moi- 
eties. Since the optical dipole of 7T7t* and of the n7r* tran- 
sition of the azobenzene moiety is directed along its long 
molecular axis, the fragments oriented perpendicular to 
the actinic light polarization vector, E, then become al- 
most inactive, whereas the other with suitable orientation 
are active undergoing photoisomerization. These trans- 
cis-trans photoisomerization cycles are accompanied by 
rotations of the azobenzene chromophores resulting even- 
tually in an orientation of the long axes of the azoben- 
zene fragments along all directions normal to the polariza- 
tion vector of the incident actinic light. Non-photoactive 
groups then undergo reorientation due to cooperative mo- 
tion or dipole interaction jn], |ll|, g |l|, Q. 

The above scenario, known as photoorientation mech- 
anism, assumes angular redistribution of the long axes of 
the trans molecules during the trans-cis-trans isomeriza- 
tion cycles. It was initially suggested in for the case, 
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when the cis state has a short lifetime, that it becomes 
temporary populated during this process but reacts imme- 
diately back to thermodynamically stable trans isomeric 
form. 

Another limiting case is known as an angular selective 
hole burning mechanism (photoselection) and occurs when 
the cis states are long living. In this case POA is caused by 
the selective depletion of the trans isomeric form during 
the establishment of the steady state Q . The anisotropy 
induced in this way is not long term stable and disappears 
as a result of the thermal back reaction. The photoorienta- 
tion process in the steady state of the photoisomerization 
takes place simultaneously, but it needs longer time to 
saturate. Generally, both mechanisms contribute to POA. 

It is evident from the foregoing that the actinic light 
results in the orientation of azobenzene chromophores per- 
pendicular to the polarization vector E. These directions 
can be thought as equivalent provided that the symmetry 
group of the system includes rotations. From the experi- 
mental results, however, the latter is not the case. In par- 
ticular, it was found that the photoinduced orientational 
structures can show biaxiality [@, [l7|, ll| . The variety 
of orientational configurations (uniaxial, biaxial, splayed) 
with different spatial orientations of the principle axes can 
be expected depending on many factors such as chemical 
structure of polymer, method of film preparation, irradi- 
ation conditions and so on. 

In the past years this spatial character of the photoori- 
entation has not received much attention. It was neglected 
in the bulk of experimental and theoretical studies of POA 
in azobenzene containing polymers [§, ||, @, ||, [| —< Bl- 
One conceivable reason for this can be the lack of appro- 
priate experimental methods. On the other hand, until re- 
cently, the problems related to the 3D orientational struc- 
tures in polymeric films has not been of major interest 
for applications. But such kind of studies are currently of 
considerable importance in the development of new com- 
pensation films for liquid crystal (LC) displays |19 and 
the pretilt angle generation by the use of photoalignment 
method of LC orientation [M ■ 

The known methods suitable for the experimental study 
of the 3D orientational distributions in polymer films can 
be divided into two groups. 

The methods of the first group are based on absorp- 
tion measurements. These methods have the indisputable 
advantage that the order parameters of various molecular 
groups can be estimated from the results of these mea- 
surements. Shortcomings of the known absorption meth- 
ods |(| [2l] are the limited field of applications and the 
strong approximations. 

The second group includes the methods dealing with 
principle refractive indices. Recently several variations of 
the prism coupling methods have been applied to measure 
the principle refractive indices in azopolymer films j2^, ^3], 
p4j . These results, however, were not used for in-depth 
analysis of such features of the spatial ordering as biaxi- 
ality and spatial orientation of the optical axes depending 
on polymer chemical structure, irradiation conditions etc. 



Our goal is a comprehensive investigation of the pecu- 
liarities of 3D orientational ordering in azopolymers. The 
present work is a part of the study focused on the orienta- 
tional biaxiality and the transition from biaxial to uniaxial 
structures caused by the polarized actinic light. 

The paper is organized as follows. 

In Sect. |^ we describe our combined approach based 
on using the methods that deal with both absorption and 
birefringence measurements. The modified null ellipsom- 
etry method is employed to study the general structure 
of the anisotropic polymer films. The components of or- 
der parameter tensor of the azobenzene chromophores are 
estimated from the results of the UV absorption measure- 
ments. 

Material of Sect. || comprises the theoretical part of the 
paper. We begin with the analysis of general kinetic rate 
equations and show how the known results [ll[] can be 
recovered by using our theoretical approach. Then we for- 
mulate the phenomenological model of the photoinduced 
ordering in azopolymers that accounts for biaxiality of the 
induced structures and long term stability of POA. After 
computing the order parameter components of azobenzene 
units for different irradiation doses we find that the pre- 
dictions of the theory are in good agreement with the data 
obtained experimentally. 

Finally in Sect. ^ we draw together the results and 
make some concluding remarks. 



2 Experimental 

2.1 Samples preparation and irradiation procedure 

We investigated POA using poly[octyl(4-hexyloxy-4'-nitro)- 
azobenzenemalonate] as model polymer which synthesis 
is described in p5[ . The thermal properties of the poly- 
mer are characterized by the transition temperatures Ci- 
32°-C 2 -44°-S 52°-N-55°-I detected by DSC and polariz- 
ing microscopy of the polymer in the bulk. (Two crys- 
talline states are labelled Ci and C2; the symbols S and N 
stand for nematic and smectic mesophases, respectively; 
the symbol I corresponds to the isotropic melt.) The poly- 
mer was solved in dichloroethane and spincoated on the 
quartz slabs. The prepared films were kept at the room 
temperature for 24 h for the evaporation of solvent. The 
thickness of the films of about 200-600 nm were measured 
with a profilometer of Tencor Instruments 

In order to induce anisotropy in the films, we used the 
irradiation of a Hg lamp in combination with an inter- 
ference filter (365 nm). The intensity of the actinic light 
was about 1.0 mW/cm 2 . A Glan-Thompson polarizer was 
applied for the polarization of the UV light. A normal in- 
cidence of the actinic light was used in our studies. 

The irradiation was provided in several steps followed 
by both birefringence end absorption measurements. In or- 
der to have ordering processes in the films completed after 
switching off the irradiation, the waiting time interval be- 
fore the measurements was longer than 15 min. According 
to |T^| , it corresponds to a time period which guarantee 
that all cis isomers react back to the trans form. 
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Fig. 1. Experimental setup for the null ellipsometry measurements. 



2.2 Null ellipsometry method 

Instead of the prism coupling methods commonly used for 
the estimation of principle refractive indices we applied 
null ellipsometry technique prj| dealing with birefringence 
components. By this means we have avoided some disad- 
vantages of prism coupling method such as the problem of 
making optical contact between the prism and the poly- 
mer layer. 

The optical scheme of our method is presented in Fig. [I . 
The polymer film is placed between crossed polarizer anc 
analyzer and a quarter wave plate with the optic axes ori- 
ented parallel to the polarization direction of the polarizer. 
The elliptically polarized beam passed through the sample 
is transformed into the linearly polarized light by means of 
the quarter wave plate. The polarization plane of this light 
is turned with respect to the polarization direction of the 
polarizer. This rotation is related to the phase retardation 
acquired by the light beam after passing through the film 
under investigation. It can be compensated by rotating 
the analyzer to the angle (f> that encodes information on 
the phase retardation. 

This method used for the normal incidence of the test- 
ing light is known as the Senarmont technique. It is suit- 
able for the in-plane birefringence measurements. 

Using oblique incidence of the testing beam we have 
extended this method for estimation of both in-plane, n y — 
n x , and out-of-plane n z — n x birefringence (n x , n y and 



n z are the principle refractive indices of the film shown 
m Fig. 0). In this case, the angle <p depends on the in- 
plane retardation (n y — n x )d, the out-of-plane retardation 
(n z — n x )d and the absolute value of a refractive index of 
the biaxial film, say, n x . 

We need to have the light coming out of the quarter 
wave plate almost linear polarized when the system ana- 
lyzes the phase shift between two orthogonal eigenmodes 
of the sample. In our experimental setup this requirement 
can be met, when the x axis, directed along the polar- 
ization vector of the actinic light, is oriented horizontally 
or vertically. Dependencies of the analyzer rotation an- 
gle 4> on the incidence angle of the testing beam were 
measured for both vertical and horizontal orientation of 
the x axis. The value of n x was measured with the Abbe 
refractometer independently. 

By using Berreman's 4x4 matrix method j27j, the 
^-dependencies of <j> were calculated. Maxwell's equations 
for the light propagation through the system of polarizer, 
sample and quarter wave plate were solved numerically 
for the different configurations of optical axes in the sam- 
ples. The measured and computed <fr versus 9 curves were 
fitted in the most probable configuration model using the 
measured value of n x . 

We conclude on alignment of the azobenzene fragments 
from the obtained values of (n y — n x )d and (n z — n x )d 
assuming that the preferred direction of these fragments 
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coincides with the direction of the largest refractive index. 
More details on the method can be found in our previous 
publication fL8f . 

In our setup designed for the null ellipsometry mea- 
surements we used a low power He-Ne laser (A = 632.8 nm), 
two Glan-Thompson polarizers mounted on rotational sta- 
ges from Oriel Corp., a quarter wave plate from Edmund 
Scientific and a sample holder mounted on the rotational 
stage. The light intensity was measured with a photodi- 
ode. The setup was automatically controlled by a personal 
computer. The rotation accuracy of the analyzer was bet- 
ter than 0.2 degree. 

2.3 Absorption measurements 

The UV/Vis absorption measurements were carried out 
using a diode array spectrometer (Polytec XDAP V2.3). 
The samples were set normally to the testing beam of a 
deuterium lamp. A Glan-Thompson prism with a computer- 
driven stepper was used for polarization of the testing 
beam. The UV spectra of the original as well as irradiated 
films were measured in the spectral range of 220-400 nm 
with the rotation step of polarizer of 5 degree. 

From these measurements the optical density compo- 
nents corresponding to the absorption maximum of azoben- 
zene chromophores were estimated for the polarization di- 
rection of the testing light parallel to x and y axes, re- 
spectively. We denote them as D x and D y , respectively. 
The out-of-plane component, D was estimated by mak- 
ing use the method proposed in || |28| . The latter assumes 
that the sample has uniaxial structure with in-plane po- 
sition of the axis of anisotropy at the instant of time to. 
It implies that D z (to) = D x (to) and the total absorption 
can be estimated as follows 

Dtotai = D x {t a ) + D y {t ) + D z {t ) = 

= 2D x (t ) + D y (to) . (1) 

When the number of trans azobenzene units does not 
change considerably, the total absorption is constant and 
the value of D z at instant of time t can be determined 
from the following equation: 

D z (t) = D total - D x (t) -D y (t), (2) 

where D x (t) and D y (t) are the experimentally measured 
parameters. 

2.4 Experimental results 

Figure ^| shows the phase shift <f> versus the incidence angle 
9 curve measured with null ellipsometry method for the 
non-irradiated polymer film. The curves (f)(8) measured for 
vertical and horizontal position of x axis overlaps. In addi- 
tion, it is seen that there is no phase shift for normal light 
incidence (6 = 0). So we arrive at the conclusion that the 
in-plane indices are matched: n y = n x . The film, however, 
possesses out of plane birefringence (n z — n x )d = —20 nm 
that results in a phase shift at oblique light incidence. 



The fitting gives the following relations for the prin- 
ciple refractive indices: n z < n x — n y and n z — n x w 
—0.1. The film, therefore, demonstrates negative birefrin- 
gence with the optical axis normal to the film surface. 
The relationship between the three indices suggests that 
the azobenzene fragments are randomly distributed in the 
plane of the film with no preferred direction for their ori- 
entation (a degenerate in-plane distribution). 

Figure ^ shows the measured (j> versus 9 curves for the 
same polymer film after 2 h of UV light irradiation. Curves 
1 and 2 correspond to vertical and horizontal position of 
the film x axis, respectively. 

According to the modeling, positive phase shift corre- 
sponds to the axis in the horizontal direction having the 
higher in-plane refractive index n y perpendicular to UV 
light polarization and the lower in-plane index n x . From 
the curve fitting we have: n y — n x = 0.25 ((n y — n x )d « 
50?im), (n z — n x )d = nm, n y > n x = n z . The light in- 
duced structure is positive uniaxial with the optical axis 
perpendicular to the UV light polarization. In this case, 
the azobenzene fragments show planar alignment perpen- 
dicular to the UV light polarization. 

The fitted values of the in-plane, n y — n Xl and out-of- 
plane birefringence, n z — n x , for various irradiation times 
are presented in Fig. |^. For small irradiation doses the 
principle refractive indices are different: n z < n x < n y . 
The in-plane birefringence monotonously increases up to 
the saturated value with the increase of irradiation time. 
On the other hand, the difference between n x and n z 
decreases and becomes negligible in the photosaturated 
state. So the film is biaxial at the intermediate stages of 
irradiation, whereas the photosaturated state can be char- 
acterized as an uniaxial structure. The optical axis of this 
structure lies in the plane of the film and is directed along 
the y axis. 

The null ellipsometry identifies general orientational 
structure, but it does not provide the means to estimate 
the order parameters of various molecular groups. Differ- 
ent absorption methods are common for this purpose. In 
this case the wavelength of testing light is tuned to the 
absorption maximum of the selected molecular fragments. 

In order to estimate the order parameter of azobenzene 
units we carried out UV absorption measurements in the 
absorption maximum of azobenzene chromophores. The 
UV/Vis spectrum of the studied azopolymer is presented 
in Fig. |] (curve 1). It contains the intensive absorption 
band with the maximum at A = 377 nm corresponding to 
the 7T7T* transition of trans azobenzene fragments. 

The spectrum reveals polarization splitting during ir- 
radiation with polarized light. The polarization compo- 
nents D x and D y , measured just after switching off the 
actinic light, are depicted in Fig. ^ as 2 and 2', respec- 
tively. These spectra show changes that become stationary 
for approximately 10 min. 

The stationary spectra D x and D y are shown in Fig. ^ 
as curves 3 and 3', respectively. In order to have the azoben- 
zene units relaxed to the stationary state, the components 
D x and D y were measured in 15 min after each irradiation 
period. 
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Fig. 2. Dependencies of the phase shift cj> on the incidence angle 9 measured before (curve 1) and after irradiation (curves 2 
and 2'). 



The experimentally measured absorption components 
D x and D y before irradiation and for different irradia- 
tion doses are presented in Fig. ||a. Kinetics of D x and 
D y is typical for reorientation mechanism of azobenzene 
units 0. Both curves reveal saturation. As it was shown 
by the null ellipsometry method, the saturated state of 
the polymer film under consideration is uniaxial with the 
in-plane orientation of the anisotropy axis. 

In order to prove that the method described in Sect. [Q| 
can be applied to estimate D z , we need to show that, 
as compared to the non-irradiated film, the number of 
trans isomers does not change in the state relaxed after 
the irradiation. This is the case under the lifetime of cis 
isomers is shorter than the time of spectral relaxation after 
switching off the actinic light. 

In order to estimate the lifetime of cis isomers we 
have measured relaxation of the spectral changes at At = 
377 nm after irradiation with non-polarized light. Inci- 
dence directions of both actinic and testing light were 
approximately normal to the film. It was found that the 
relaxation curve D(t) contains two components with char- 
acteristic times of 1.0 sec and 4 min, respectively. The first 
value can be attributed to cis-trans transition of azoben- 
zene chromophores. The other time can be related to ei- 
ther the small fraction of long living cis isomers or more 
likely to the orientational relaxation of the azobenzene 



units. In any case, the waiting time of 15 min ensures that 
we definitely have all the cis isomers transformed into the 
trans form. 

So in our measurements the concentration of trans iso- 
mers is preserved and the method described in Sect. 2J? 
can be applied. The values of D z calculated for various ir- 
radiation doses using Eqs. ([!])- (|^) are presented in Fig. 
Dependencies D x (t), D y {t) and D z (t) show that the pho- 
toinduced ordering is mainly due to the in-plane reorien- 
tation of azobenzene fragments in the y direction. In ad- 
dition, slight reorientation in the z direction is observed. 

The orientational structure is generally described by 
the tensor Sy , which is diagonal when the coordinate axes 
directed along the principle axes of the film. The diagonal 
elements S xx = S x , S yy = S y and S zz = S x are related 
'Sorption cor 

example, 



to the absorption components D x , D y and D 



For 



2D X - (D y + D z 
2{D X +D y + D z 



(3) 



The components S y and S z can be obtained by the cyclic 
permutation in the expression (^|). 

The values of S x , S y and S z calculated using equa- 
tion ([|) are presented in Fig. |5|b. As is seen, during the 
initial stage of irradiation orientational configuration is bi- 
axial, whereas the initial and the photosaturated states are 
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Fig. 3. In-plane and out-of-plane birefringence as a function of the irradiation time. 



uniaxial. These structures are characterized by the follow- 
ing order parameters: S x = S y « 0.035 and S z « —0.07 
for the non-irradiated film; S x — S z ~ —0.1255 and S y ss 
0.251 for the film in the photosaturated state. 

The transition between biaxial and uniaxial photoin- 
duced orientations will be the subject of subsequent stud- 
ies. We believe that a tendency to form an uniaxial struc- 
ture is related to an intrinsic property of self-organization 
of mesogenic groups. Theoretically, similar tendency of ne- 
matic liquid crystals can be related to some specific fea- 
tures of the phase transition reflected in the form of the 
mean field free energy |29|. In this connection we assume 
that the irradiation with polarized light causes both the 
photoorientation of azochromophores and the self-aggre- 
gation as it was found for other liquid crystalline azopoly- 
mers |50l |3li. 



3 Theory 

Before we proceed with theoretical considerations let us 
emphasize the following distinguishing features of POA in 
liquid crystalline azopolymers: 

(a) in contrast with the reversible POA, considered in jl5], 
|32| , H||, HJ, p5|] , where POA disappears after switching 
off the irradiation, POA in liquid crystalline azopoly- 
mers can only be thermally erased by heating above 



clearing temperature. The long term stability of POA 
is caused by the photoorientation of the azobenzene; 
(b) the biaxiality effects discussed in Sect. 0. 

Clearly, we can conclude that the photoorientation is a 
non-equilibrium process in a rather complex polymer sys- 
tem and it still remains a challenge to develop a tractable 
theory treating all the above points adequately. 

As far as the long term stability of POA is concerned, 
the reorientation of the azobenzene groups can be assumed 
to result in the appearance of a self-consistent anisotropic 
field that support induced anisotropy. This field comes 
from anisotropic interactions between the azobenzene frag- 
ments and rearrangement of the main chains and other 
non-absorbing fragments. In other words, the photoin- 
duced orientational structures can be regarded as a result 
of the photo-reorientation and thermotropic self-organiza- 
tion processes |kJ. 

There are two phenomenological models based on sim- 
ilar assumptions: the multidomain model proposed in H 
and the model with additional order parameter, attributed 
to the polymer backbone and introduced to make the 
steady state degenerate fuJl . 

Despite these models look different it is clear that 
they incorporate the long term stability by introducing 
additional degree of freedom (subsystem) which kinetics 
would reflect cooperative motion and account for non- 
equilibrium behavior. 
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Fig. 4. UV absorption spectra measured before (curve 1) and immediately after the irradiation over 1 h (curves 2 and 2'). The 
stationary spectra D x and D y (the waiting time is 15 min) are shown as the curves 3 and 3', respectively. 



In this section we describe our theoretical approach to 
the kinetics of the photoinduced reorientation. We begin 
with the analysis of general master equations and spe- 
cialize then the rates of the involved transitions. The re- 
sulting kinetic equations for order parameters are derived 
after making assumptions on the form of angular redistri- 
bution probabilities and the order parameter correlation 
functions. In addition, we show that this approach can 
be employed to derive the known phenomenological mod- 
els H, 0], ^2) . Finally at the end of the section we present 
the numerical results obtained by solving the kinetic equa- 
tions for the order parameters and concentrations. 



3.1 Master equations 

We shall assume that the dye molecules in the ground state 
are of trans form with the orientation of the molecular axis 
defined by the unit vector n. The latter is specified by the 
polar, 9, and azimuthal, (p, angles: n =(sin 6 cos 0,sin 9 sin <f>. 
cos 8). 

Angular distribution of the trans molecules at time 
t is characterized by the number distribution function 
N tr (n,t). Molecules in the excited state have the cis con- 
formation and the corresponding function is N C i S (n,t). 



Then for the number of trans and cis molecules we have 



N tr (t) = Nn tr (t) = / i\T tr (n,f)dn, 



(4) 



N ci8 (t) = Nn cls (t) = J N cts {h 7 t)dn 7 (5) 
N = N tr (t) + Nd a (t) , n tr (t)+n cis (t) = l, (6) 



where 



dn 



sin^d^ and N is the total 



number of molecules. 

We shall refer the additional subsystem that is able to 
accumulate induced ordering of the side chain molecules 
as a polymer system (matrix). From the phenomenologi- 
cal point of view, this system can be thought to represent 
some collective degrees of freedom of non-absorbing units 
such as main chains. We shall suppose that it is character- 
ized by the angular distribution function / p (n,i), so that 
N p (n, t) = A r p /p(n, t). Note that the coefficient N p can be 
considered as an effective number of the units related to 
the polymer system. But, more precisely, this factor de- 
termine the relations between different thermal relaxation 
constants (see Eq. (jljj)). 

The starting point of our approach is the kinetic rate 
equations taken in the following general form |36|, [37j : 



8 



Yu. Zakrevskyy, O. Yaroshchuk, A.D. Kiselev et al: Spatial reorientation in LC azopolymer 



0.9 



0.8 



S 
d 

y 
S 



0.7 



°0.6j 



t>2 



0.5 



0.4 



JET 



■ 



JET 



1 


1 1 1 








□ D y (exp.) 


G — 


O D x (exp.) 


0-~ 


O D z (exp.) 



(a) 



50 100 
Time, min 



150 




,,-Er"" 


-n — — 3 




— theory 


G — 


-O S x (exp.) 


□ — 


-□ S y (exp.) 




-O S z (exp.) 


(b) 






-@— © 



0.3 



0.2 



2 

o.i! 

On 
| 



-0.1 



50 100 150 
Time, min 



200 



0-2 
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time. Theoretical curves for the diagonal components of the order parameter tensor are shown as solid lines. 



dN a 
dt 



dt 



Diff /3#Q 



]T J [w a ^(n,h')N (h\t) 



W^ a (a',n)N a (a,t) 



da' , 



(7) 



where a ,0 S {tr, cis, p}. 

The first term on the right hand side of Eq. (0) is due to 
rotational diffusion of molecules in trans (a = tr) and cis 
(a = cis) conformations. Note that the terms proportional 
to W a ^ a can be incorporated into this diffusion term. In 
what follows it is supposed that 



~dN tr ' 




fd/p" 




dt 


Diff 


dt 


Diff 



= 0. 



(8) 



Now in order to proceed we need to specify the rates 
of the transitions. 



3.2 Transition rates 



The trans-cis transition is stimulated by the incident UV- 
light quasiresonant to the transition. Assuming that the 
electromagnetic wave is linearly polarized along the x— 



axis, the transition rate can be written as follows p3| , p5[ : 

W cis ^ tr (a,a') = r c _ t (n,n')Ptr(fi')> 

y j r c _ f (n,n')dn= 1, (9) 

p tr (ti) = M-^r^xXf^')^; - 

= q t I{\ + un 2 x ) = q t I{\ + u {2S X + l)/3) , (10) 

where cr^ tr \a) is the tensor of absorption cross section 
for the trans molecule oriented along ft: afp = er^ <% + 

a { t r) )/o [ t r) is the ab 



a 



(tr) 



) rii n 



j ■ 



- I (tr) 

u = (a\ { 



sorption anisotropy parameter; Hu)t is the photon energy; 
$tr^cis is the quantum yield of the process and r tr (ft, ft') 
describes the angular redistribution of the molecules ex- 
cited in the cis state; / is the pumping intensity. 

Similar line of reasoning applies to the cis-trans tran- 
sition, so we have 



w tr ^ cls (n, n') = lc r£ p ] (ft, ft') + q j r t ( Z a > (ft, ft'), 

q c = $cis^trans<T {ciS) , (H) 

where 7 C = l/r c , r c is the lifetime of cis molecule and the 
anisotropic part of the absorption cross section is disre- 



garded, a\^ s) = cr^ ls) = . Eq. (0) implies that the 
process of angular redistribution for induced and sponta- 
neous transitions can differ. Note that the normalization 
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condition for all the angular redistribution probability in- 
tensities is 



J r Q _ /3 (n,n')dn= 1 



(12) 



The remaining part of transitions describes equilibrat- 
ing between the side chain absorbing molecules and the 
polymer system. The corresponding rates can be taken in 
the form: 



W a ^ p (n, n ) = 7 Q _ p r Q _ p (n, n ), 

Wp^ a (n,n') = 7 J) _ a r j) _ a (n,n'), a £ {tr, cis} , (13) 



where 7 Q - P and ja-p are angular independent. In addi- 
tion, since thermal relaxation does not change the number 
of molecules in a particular state we have the relation for 

7 a _ J) and 7 p _ a : 



iV p 7 Q 



JVn Q 7 P „ Q 



(14) 



As it was mentioned above, this equation relates the ther- 
mal relaxation constants of the polymer and the frag ments 
through the coefficient N p , introduced in Sect. 3.1. 



3.3 Model 



At this stage it is convenient to introduce normalized an- 
gular distribution functions, f a (n,t): 



N a (h,t)=Nn a (t)f a (h,t). 



(15) 



From Eqs. (0), (g) and ( |Tl| ) it is not difficult to obtain 
equation for nt r (t): 



dn tr (t) 
dt 



= (7c + qJ) n cls (t) - (Ptr)trntr(t) , (16) 



where the angular brackets (. . . ) a stand for averaging over 
angles with the distribution function f a . Owing to the 
condition (|12j), this equation does not depend on the form 
of the angular redistribution probabilities. 

From the results of the previous section and Eq. ( |l6|) 
we derive the equations for the distribution functions 



n C i s (t) — = n cls (t)< 



dfc 



dt 



Diff 



7cis/c lS (n, t) ) - (P t r}trn tr (t) + 



+ n tr (t) J r c - t (h,n')P tr (n')f tr (h',t) dn' + 

p(n, n')/ p (n', t) dn' , (17) 



n tr {t) — g"' - = -n tr (t) [Ptr(n) - (Ptr)tr + 7tr] /tr(n, t) 
+ 7c n cls (t) J r t ^ (n, n')/ cis (n', t) dn' + 

+ qjn cis (t) J r t {in c d) (n, n')f cis (n', t) dn' - 

- (7c + q c P) n cis (t) / tr (n, t) + 

+ lt r n tr {t) J r t _ p (n,n')/ p (n',i)dn', (18) 

- J r p _ t (fl, n')/ tr (n', *) dfi'J - 7p " s) n cis (t) ^/ p (n, t) - 



y r p _ c (n ) n / )/ci s (n / ,f)dn / J , 



(19) 



where j tr = Jp-trans, Ids = Ip-as and 7^ = Nj a /N p . 
These equations supplemented with Eq. ( J16| ) are derived 
on the basis of quite general considerations. They can be 
regarded as a starting point for the formulation of a num- 
ber of phenomenological models. We can now describe our 
model. 

Our basic assumptions of the angular redistribution 
operators r a ~i3 are as follows 

lcis = 0, r£ p } (n, n') = r t _„(fi, ft') = 

= r p _ t (fi,fi') =5(n-fi'), (20) 
r t ( _" c d) (n,n') = f tr (n,t), r c _ t (fi,fi') = f cls (n,t) . (21) 

It gives the resulting system of kinetic equations: 



dfc 



dt 



dfc 



dt 



Dill' 



Hi 



dt 

+ 7c rids ( /, 



{{Ptr)tr - Ptr - 7tr) n tr ftr + 
far) + 7tr «tr /p , 



^ 7p n tr(ftr fp) 
Atr) 



(22a) 

(22b) 
(22c) 



where 7 P = 7 P 

Clearly, the meaning of Eq. (|20|) is that the molecules 
do not change orientation under spontaneous transitions. 
On the other hand, from Eq. ( pl| ) projecting onto the 
angular distribution function of the corresponding state 
describes angular redistribution for the stimulated transi- 
tions. 

In order to explain the meaning of the projectors, note 
that the multidomain model considered in |1 can be de- 
rived from Eqs. (Q) - (Q by putting 7 C = 0, 7 CIS = j tr 
and assuming that all angular redistribution probabilities 
are equal to the equilibrium distribution, r a -/3 (n, n') = 
p(n), determined by the mean field potential W(h): p(h) oc 
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(c) t = 120.0 min 



Fig. 6. The orientational distribution functions of trans molecules at different illumination doses. These surfaces are defined 
by the equation: r(9.cf>) = 4tv nt r (t) ftr(d, 4>, i) in the spherical coordinate system at irradiation time t. 



exp(— W). In other words, this procedure introduces the 
mean field potential by assuming that the angular redistri- 
bution operators r a -p act as projectors onto the equilib- 
rium distribution. Note that the results of JO, [l2) cor- 
respond to the case where r^ p J (n, n') = /tr(n, t) and 
ri!" c d) (n,n') = (47r)- 1 . 



3.4 Order parameters 

We can now deduce the equations that describe the tem- 
poral evolution of the diagonal components of the order 
parameter tensor 



S(n) 



(23) 



The components of interest can be expressed in terms 
of Wigner /^-functions |58| as follows 

S x = (3n 2 x - l)/2 = 

= -(^ (n)-V6Re^ (n))/2, 



(24) 



S y = (3n 2 y - l)/2 = 

-(D 2 Q0 (n) + V6RcD 2 20 (h))/2, 

S z = (3n 2 z -l)/2 = D 2 oa (h). 



(25) 
(26) 

The simplest ca se o ccurs for the order parameters of 
cis molecules. Eq. (22a) yields the following result: 



no( cis ) 

ij _ apt q( cis ) 



(27) 



Yu. Zakrevskyy, O. Yaroshchuk, A.D. Kiselei 



where = (Sij(n)) a and D r is the rotational diffusion 
constant. Clearly, our assumptions correspond to the case 
where the presence of cis molecules is of minor importance 
for ordering kinetics. 
Eqs. 



fllPf), ( f22b| ) and (22c) give the following system 
for the components of the order parameter tensor: 



OS. 



itr) 



THr 



-2/3 q t Iu n tr G 



(tr) 
ij\ xx 



+ l c n cls (st ] ' S^ r) ) + lt r n tr (S^ - S^) 



dS 



(p) 



at 



- lp n tr (S^ - S^) . 



(28a) 
(28b) 



where G^ mn is the correlation function of the order pa- 
rameter components S# (n) and S mn (n) that is defined by 
the following relation 



(29) 



Computing the order parameter correlation functions 
that enter Eqs. ( |28| ) requires the knowledge of details on 
microscopic interactions and, in general, for nonequilib- 
rium system it can be rather involved and sophisticated. 
In this paper, we shall adopt the simplest "kinematic" pro- 
cedure to approximate the correlators. It implies that after 
writing the products of D-functions as a sum of spheri- 
cal harmonics we neglect the high order harmonics with 
angular momentum j > 2. In particular, we have 



(S2)fr«l/5 + 2/7(S a ) tri ae{x, y, z} . 



(30) 



Applying this procedure to Eqs. fl2q) leads to the result 
given by 



n tr — = 2u/3q t I(l/5 + 2/7 S-S 2 )n tr 

- IcflcisS + +^t r ntr(Sp - S) , 



OAS 

~~dT 



-2u I '3 q t I "(2/7 + S)n tr AS - 



j c n cls AS + ~f tr n tr {AS p - AS) 



dSp 
dt 



dASp 
dt 



-l P n tr (S p - S) , 



-j p n tr (AS p - AS) , 



(32) 



(33) 



(34) 



where S = (S y ) t r, AS = (S x - S x ) tr , S p = (S y ) p and 
AS P = (S x -S z ) p . Eqs. @ - @ combined with Eq. @ 
form the system of kinetic equations for our phenomcno- 
logical model. 
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Fig. 7. Initial orientational distribution at section by the x-y 
plane is isotropic. It assumes anisotropy along the y axis as 
illumination time increases. 



3.5 Numerical results 



Theoretical curves depicted in Fig. [|b are calculated by 
solving the equations deduced in the previous section. This 
procedure involves computing the stationary values of S 
and AS to which the order parameters decays after switch- 
ing off the irradiation at time t. In addition, we need to 
take into consideration the difference between the order 
paramet ers defined by Eq. (|3|) and the order parameters 
of Sect. Since A oc (1 + tt(2Sj + l)/3), these order 
parameters differ by the factor it/(3 + u). 

According to the experimental data, the lifetime of cis 
molecules r c (<y c = l/r c ) is about 1.0 sec, whereas the re- 
laxation time after switching off the irradiation can be 
estimated at 4 min. Since the theoretical value of this 
relaxation time is 1/(7 P + Jtr)i the relaxation times r p 
(7p = 1/tj,) and r tr (jtr = 1/ T tr) can be taken to be equal 
min. 

We estimated the absorption cross sections <j( cls ) and 
<j( tr ) f rom the UV spectra of the polymer solved in toluene. 
These spectra were measured before and during irradi- 
ation. In the latter case the solution was in the photo- 
saturated state. The absorption bands were then decom- 
posed into the bands of trans and cis isomers to yield 
the corresponding values of the extinction coefficients at 
A f = 365 nm. In order to compute these coefficients we 
followed the procedure described in Q. The resulting 
estimates can be written as follows: (fiw t ) _1 er( cls V c 7 sa 
0.2 x 10" 2 and a^/a^ « 2.5 for I = 1 mW/cm 2 , 



where cr ( *^ = (aj,* r) + 2<r [ "'>)/3 is the average absorption 
cross section of the trans fragments. 

Then given the quantum yield <P C i S ^trans and the ex- 
perimental value of the order parameter in the photo- 
steady state, S st = 0.251 (1 + 3/u), we can compute the 
anisotropy parameter u from the equation for S s t- The 
latter can be derived from Eqs. (0) and ([n|) by setting 
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the time derivatives on the left hand sides equal to zero: 

7c5 rt (l + «(l-5 rt )/3) = 

= 2u/3 ( 7c + q c I) (1/5 + 2/7 S a t - S*) . (35) 

The theoretical curves in Fig. ||b are calculated at 
+trans = 10% and ^ds^trans = 5% that is to yield the 



8.9. Note that the quan- 



value of the ratio a 
turn efficiencies are of the same order of magnitude as the 
experimental values for other azobenzene compounds fio] ] . 
On the other hand, we have q c w 0.2 x 10 -3 cm 2 /mJ that 
is about the value given in ]^]. 

The computed order parameter components can be 
used to illustrate orientational distributions of trans frag- 
ments maintained after different illumination doses. Fig. |6| 
shows the surfaces r(6.<f>) = 4ir nt r (t) ftr(6, t) that indi- 
cate the angular redistribution in the course of irradiation. 
Note that we have truncated the expansion for the distri- 
bution function f tr by neglecting the high order spherical 
harmonics: 



47r/ tr (fM) ~1 



+ 



-5 

(S x )tr 



(S z ) tr D* (n) 



-^k 2 ReA 2 o( 



(36) 



Sections of the surfaces depicted in Fig. || by the x- 
y and x-z coordinate planes are shown in Figs. ^ - 0. 
As is seen from Fig. [?], the angular distribution in the 
x-y plane becomes anisotropic under the action of light, 
whereas Fig. || indicates that the distribution in the x-z 
plane goes isotropic. 



tsj o - 



t = 0.0 min 

t = 10.0 min 

t = 120 min 




o 

X 



Fig. 8. Increase in the illumination dose renders the orienta- 
tional distribution of the trans fragments at section by the x-z 
plane isotropic. 



4 Conclusions 

In this paper we demonstrate that the combination of the 
UV/Vis absorption spectroscopy and the suitably modi- 
fied null ellipsometry method is a tool appropriate for a 



comprehensive study of POA in films of liquid crystalline 
polymers with azobenzene side groups. 

We found that initially the spincoated films under in- 
vestigation are characterized by a preferred in-plane orien- 
tation of the azobenzene fragments, whereas the long axes 
are randomly distributed in this plane. So, the resulting 
structure is isotropic in this plane, but in the space it is 
uniaxial with the optical axis normal to the film surface 
resulting in a negative birefringence (see Sect. |2.4| ). 

Increasing irradiation doses results in an anisotropic 
order maintained in the film after switching off the lin- 
early polarized UV light. This photo- induced structure 
corresponds to a biaxial orientational order of the azoben- 
zene moieties. But it becomes uniaxial at sufficiently large 
doses reaching the photosaturated state of the photoori- 
entation process. The process results usually in an oblate 
order in the case of amorphous polymers. The uniaxiality 
in the studied case could be related to liquid crystalline 
character of the polymer or the initial order of the film 
which act in combination with the linearly polarized light 
as aligning force resulting in an uniaxial in-plane order. So, 
the action of the actinic light can be regarded as a factor 
stimulating both photoorientation and thermotropic self- 
organization. 

Quantitatively, experimental results on the kinetics of 
the photoreorientation were descri bed in terms of the or- 
der parameters introduced in Sect. 2.4. It makes the com- 
parison between experimental and theoretical results rel- 
atively direct. On the other hand, it raises the question as 
to the validity of the procedure used to determine the out- 
of-plane absorption D z . In our case this can be justified 
for the number of trans fragments is shown to remain un- 
changed in all the relaxed states caused by the lifetime of 
cis isomer. So, the conclusion is that the photoreorienta- 
tion in the film occurs through the mechanism of angular 
redistribution. Note that, as a matter of fact, the latter 
was implicitly as sumed in our theoretical model formu- 
lated in Sect. B.3. 



Despite the theoretical considerations of Sect. |3| are 
rather phenomenological they emphasize the key points 
that should be addressed by such kind of theorie s. T hese 
are the angular redistribution probabilities (Sect . |3.3[) and 
the order parameter correlation functions (Sect. p.4j ). The 
redistribution operators, in particular, define how the sys- 
tem relaxes after switching off the irradiation and can 
serve to introduce self-consistent fields. The correlators, 
roughly speaking, mainly determine the character of pho- 
toreorientation and the properties of the photosaturated 
state. 

The assumptions taken in this paper give the model 
with the kinetics governed by the mechanism of angular 
redistribution. The form of the approximate correlators 
influences it in such a way that the out-of-plane reorien- 
tation is appeared to be suppressed. 

Our simple model depends on a few parameters that 
enter the equations and that can be estimated from the 
experimental data. Only the anisotropy parameter and the 
quantum yields are derived by making comparison be- 
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tween the experimental data and the theoretical depen- 
dencies. 

This theory predicts that the fraction of cis isomers 
is negligible and the orientational kinetics of the cis frag- 
ments governed by Eq. ( |27j ) is irrelevant. The latter is why 
the liquid crystalline ordering effects |4l| that would com- 
plicate the kinetics of the cis molecules can be safely omit- 
ted in Eq. (^). Note, however, that it is rather straight- 
forward to modify the model for the case where the cis 
fragments would affect the kinetics considerably We shall 
extend on the subject elsewhere. 

More detailed microscopic theoretical approach is ap- 
parently beyond the scope of the paper. Theoretically, it 
would be interesting to put this theory into the context 
of ergodicity breaking transitions [ [f2[ . This work is under 
progress. 

We acknowledge financial support from CRDF under the grant 
UP1-2121B. We also thank Dr. T. Sergan and Prof. J. Kelly 
from Kent State University for assistance with processing the 
data of the null ellipsometry measurements and helpful discus- 
sions. 
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